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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 2878 


COMBINED EFFECT OF DAMPING SCREENS AND STREAM 
CONVERGENCE ON TURBULENCE 


By Maurice Tucker 


SUMMARY 


An analysis is presented of the combined effect of a series of 
damping screens followed by an axisymmetric-stream convergence (or diver- 
gence) upon the mean-square fluctuation-velocity intensities, scales, 
correlations, and one-dimensional spectra of a turbulence field con- 
vected by a main stream. The treatment is restricted to negligible tur- 
bulence decay and linearized by postulating small fluctuation velocities 
and velocity gradients, and absence of viscosity except as simulated by 
the idealized screen action. Compressibility of the main stream is 
allowed for during passage through the contracting section. The density 
fluctuations associated with the turbulence field are regarded as 
negligible. 


Numerical results for the statistical quantities describing the 
turbulence field downstream of a screen-contraction configuration are 
obtained for the case of upstream isotropic turbulence. The action of 
the damping screens and the stream convergence is to distort this ini- 
tially isotropic field into a field of turbulence symmetric about the 
longitudinal direction with the lateral fluctuation velocities greater 
in magnitude than the longitudinal velocities. 


An approximate method of taking into account the effects of tur- 
bulence decay upon the mean-square fluctuation velocities obtained for 
the case of negligible decay is presented. This method of correction 
together with the tabulation of fluctuation-velocity ratios over an 
extensive range of conditions should prove useful for engineering 
applications. 


INTRODUCTION 
The use of fine-mesh or damping screens located in a low-speed 


settling chamber followed by a contracting passage (entrance cone) to 
attain a low-turbulence test-section flow is well known from the 
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qualitative standpoint. Dryden and Schubauer (reference 1) have pre- 
sented experimental data regarding the combined effect of screens and a 
contraction on the intensity of turbulence. . Existing theoretical studies 
are confined to either the effect of the screens or of the stream con- 
traction on turbulence. Taylor and Batchelor (reference 2) have obtained 
the effect of a damping screen located in a constant-area passage upon 

& triple Fourier integral representation of a turbulent field. The 
effect of a contraction upon a similar representation is analyzed in 
reference 3. In both references 2 and 3 initial isotropy is postulated 
in order to obtain numerical results. 


The analyses of references 2 and 3 indicate that in the absence of 
decay effects (dissipation and mixing) an initially isotropic turbulence 
field will be distorted into a field of turbulence axisymmetric about 
the mean flow direction upon passage through either a damping screen or 
an axisymmetric contraction (contraction with all cross sections simi- 
lar). An analysis of axisymmetric turbulence is given in reference 4. 
In conventional wind-tunnel configurations, turbulence that is initially 
isotropic will thus have been distorted into axisymmetric turbulence 
after passage through the first of the several damping screens and will 
remain axisymmetric while traversing the remaining screens and the fol- 
lowing contraction. Inasmuch as the expressions obtained in reference 3 
for the downstream mean-square velocity fluctuations require that the 
turbulence upstream of the contraction be isotropic, the results of ref- 
erences 2 and 3 cannot be combined in any simple manner to obtain the 
joint effect of screens and a contraction on turbulence that is ini- 


tially isotropic. 


The present analysis treats the combined effect of a series of N 
(symbols are defined in appendix A) identical damping screens and a 
downstream axisymmetric contraction upon the longitudinal and lateral 
turbulence velocity fluctuations, scales, correlations, and spectra of 
a turbulence field described by a triple Fourier integral. The config- 
uration is shown schematically in figure 1. Although compressibility 
of the main stream is allowed for during passage through the contraction, 
the density fluctuations associated with the turbulence are regarded as 
negligible. The assumption of small turbulent velocity fluctuations and 
‘velocity gradients together with the postulated absence of viscosity, as 
in references 2 and 3, implies the absence of turbulent decay processes 
and linearizes the governing equations for both the screen and con- 
traction effects. 


After a discussion of the spectrum concepts used in the present 
analysis, the preliminary portions of the analysis which borrow from the 
results of references 2 and 3 are concerned with the effect of & screen 
and of a stream contraction upon a representative wave or Fourier com- 
ponent. Briefly, the screen affects only the amplitude vector of the 
wave; the contraction acts to change both the amplitude and wave-number 
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vectors. In view of the linearized analysis and the resulting absence 
of modulation or mutual interference between the array of plane waves 
making up a field of turbulence, the correlation tensor is developed 
from the results obtained for a typical wave. The spectral tensor is 
Obtained as the Fourier transform of the correlation tensor. Turbulence 
velocity and scale ratios obtained from the spectral densities (diagonal 
components of the spectral tensor) are then given in tabular form for 
the condition of upstream isotropic turbulence. The one-dimensional 
spectra and the correlation-coefficient curves for a special case of 
upstream isotropic turbulence are also determined. Ап approximation for 
taking into account decay effects is suggested. This investigation was 
conducted at the NACA Lewis laboratory. 


ANALYSIS FOR NEGLIGIBLE DECAY 
Spectral Representation of Turbulence 


Turbulence is often regarded as an assembly of eddies of randomly 
varying size and intensity. The present analysis treats the turbulent 
Pield as a spectrum of plane sinusoidal waves with all possible wave- 
lengths, wave-front orientations, and phases. This superposition provides 
the necessary three-dimensional character to the turbulence representa- 
tion. Large eddies thus are represented by waves of large wavelength 
(small wave number). The fluctuation-velocity components a (v = 1,2, 2) 


are represented at a given instant by the triple Fourier integral 


q (x) = - [If aye — 7 dk, dk dl (1) 


where x is a position vector, Qy a wave-amplitude vector (reference 3), 
and К a Wave-number vector normal to the wave front. In order that 

the wave amplitude vector Qy be finite, the field of turbulence 
described by equation (1) is assumed to occupy a bounded region and to 
vanish everywhere outside this region. For the case treated herein in 
which the fluctuation components are related by the incompressible-flow 
form of the continuity equation 


Let | (2) 


the plane waves of equation (1) are transverse. TIn the summation of 
equation (2) the index y covers the range of values 1, 2, 3. 
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In order to obtain the spectral tensor and, in turn, the mean- 
square velocity fluctuations, it will be convenient to discuss first the 
correlation tensor and indicate its relation with the spectral tensor. 
The correlation tensor Ry s5(r) is defined as the spatial mean value of 
the коне: of the velocity component q- at x and the velocity com- 
ponent g4! at x'=x+r as x varies and the separation vector г 
of the two points remains fixed during the averaging. If it is assumed 
that the fieid of turbulence is homogeneous and statistically steady 
and that the field is confined to a parallelepiped of edges 211, cDo, 
2Ds and vanishes everywhere outside, the space average is derived in 
reference 3 as 


Ars = з "WE — Q (Q (e = ш лаш 


where Qs * (E) is the €" conjugate of Qs(k) and T is the volume 
ax ж 

80-0205 of the parallelepiped. The expression lim T Q (E)Qs (К) 

is equivalent to the spectral tensor Ds (E) defined in reference 5 as 

the Fourier transform of the correlation tensor R,»(r) 


-ik- 
BEyg(r) = ШЕ án dk, dk, dks 


| | 
ГБ (® = im — 9, (2909 Q9 (3) 


A knowledge of the spectral tensor permits, as will be shown, determi- 
nation of the various statistical quantities describing a turbulence 
field. Equation (3), which relates the spectral tensor to the wave- 
amplitude vector obtained for a typical Fourier component in the absence 
of any modulation effects,is thus basic to the present analysis. 


or 


For isotropic homogeneous turbulence fields wherein the incompress- 
ible flow form of the continuity equation is satisfied, Batchelor (ref- 
erence 5) has shown that the spectral tensor can be written 


TQ = ex) (es 2:3 n 


ne КР = ky? + kp + kg“, 5,5 = 1 for y= 5, and бүр = О for 
Y #68. 
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In matrix form 
K2 kg. -k kp aes 
ГЮ = G(k) | dk ёғ к, (4) 
Mis E E 


It is clear from the definition of the correlation tensor that for 
т = О the diagonal elements of the tensor yield the mean-square velocity 
fluctuations. In terms of the corresponding elements of the spectral 
tensor (energy spectral densities) 
co 


42 = J Ili Dr (1%) акуакак„ (s) 


The mean-square velocity fluctuations of equation (5) refer to 
spatial averages. Hot-wire instrumentation used to obtain these fluc- 
tuations, however, provides only time averages. Taylor (reference 6) 
was able to show that the spectrum of the velocity fluctuations in time 
is the Fourier transform of the spatial correlation function. Taylor's 
hypothesis (reference 7) that the main stream carries along the pattern 
of a weak field of turbulence unchanged past the point of measurement 
permits analysis of the hot-wire output signal in the form of a one- 
dimensional spectrum defined in the equivalent of spatial terms. The 
relation between the one-dimensional spectral densities Er and the 
three-dimensional spectral densities Түү is easily shown by writing 
equation (5) as 


gê = f all Dry (E) акрак, |аку = рва аку (6) 


The various statistical quantities which characterize a field of 
turbulence may be obtained from the one-dimensional spectral densities 


ав’ discussed in reference 8. Noting that == non dk, = 1, 
O 


the correlation coefficients are given by 


r4,0,0 c 
Ry (1) = Ча | Fy(K1) cos Кугу аку (7) 


E" 


Ss — —— m — À———— n. TF ————————Ó—M———————— —á————— "——-— T 
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The Fourier transform relations yield 


` 


29,2 E 
Ej (k) к. o | Ry cos Кугу dry (8) 


Two sets of characteristic lengths are customarily defined for a tur- 
bulence field. The turbulence microscales A. (mean lengths weighted 


in favor of the small eddies which are responsible for the greater part 
of the viscous dissipation) are given by 


IR 5 
3, =-(— 5) . | куёв (ку) аку (9) 
^v l/rj-0 Ф Jo 


The turbulence scales Ly (mean lengths representative of the average 
size of all the eddies) are obtained as 


L. = ne dr, = = [509]. , (10) 


day 


This physical meaning for the scale of turbulence is only applicable 
when Ry, > 0 as гу >0- 


Plane-Wave Analysis for Damping Screens 


The preceding equations indicate that the statistical quantities 
describing a field of turbulence may be obtained from the spectral 
tensor of equation (3), which is presented in terms of the plane-wave 
amplitude vectors Qy(k). The assumptions of small turbulent velocity 
fluctuations and of inviscid flow, with regard to both the main stream 
and the turbulence field convected by the main stream, linearize the 
equations which govern the action of the screens and of the contraction. 
In the resulting absence of any modulation or interaction effects 
between waves, the analysis is simplified by first treating the effect 
of a screen and a stream convergence (or divergence) upon a represent- 
ative plane wave. Superposition is then used to obtain the combination 
of these effects upon the complete assembly of plane waves which 
describes the turbulent field. 


The action of a fine-mesh or damping screen on a disturbance con- 
vected by a low-speed uniform stream may be characterized by two param- 
eters K and а. The parameter K is defined in terms of the pres- 
sure drop AP required to drive fluid of density р and velocity U 
through the screen 
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IJ pu? 


The parameter a which takes into account the side force per unit area 
was introduced by Taylor in reference 9 and relates the angles of flow 
incidence үу and flow emergence Yo shown in figure 2. It has beén 


shown experimentally that the ratio tan yo/tan y, tends to a finite 


limit o as Y], which is usually very small, tends toward zero. For 
incompressible flow the continuity equation requires that the longitu- 
dinal velocity component be unchanged after passage through the screen. 
From kinematical considerations, at the screen the ratio of downstream 
to upstream lateral velocity components equals a for small values of 
the flow incidence angle у] - 


АБ in reference 2 the uniform stream is regarded as incompressible 
and inviscid throughout the constant-area settling chamber in which the 
screens are located (station A to station B of fig. 1). A screen will, 
in general, decrease turbulent motions of larger scale than the mesh 
size and introduce turbulence of smaller scale. In the analysis the 
damping screens are assumed not to generate any wake turbulence, which 
implies that the screen mesh size and wire diameter are very small rel- 
ative to the scale of the upstream turbulence. Far upstream of the 
screen, at station A, a single plane wave carried along by the main 
stream of velocity U in the xj-direction will be designated 


~A qA ei(k- x-k,Ut) 


oe 


Coordinate axes are fixed, with the origin located at the screen and the 
positive xj-axis pointing downstream. It is shown in reference 2 on the 
basis of a steady-state disturbance analysis that far downstream of the 
Screen, at station B, the wave is transformed to 


ав av. US Ue) 


In order to satisfy conditions at the screen, it is necessary to postu- 
late disturbance fields upstream and downstream of the screen which are 
induced by the screen. These disturbance fields attenuate, vanishing 
at stations A and B. Taylor and Batchelor represent these induced 
velocities in terms of potential flows. With the velocity components 
u, v, w Of figure 2 designating the combined effect of the turbulent 
velocity fluctuations and the induced velocities, the following condi- 
tions are imposed at the screen (x4 = 0) 


Demo = (Ш - 


“< — a —  —  .— —— — O.A ` р ——<- — aea a 
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B A 
(v,w) ху=0 = a(v,w) ху=0 


The root-mean-square fluctuation velocities are taken to be small rela- 
tive to the stream velocity so that the equations of motion can be 
linearized. A further condition is imposed that the local pressure drop 
across the screen is determined by the local longitudinal velocity and the 
screen pressure-drop coefficient K. The basic relations describing 

this idealized action of a damping screen on a representative plane wave 
are then obtained in reference 2 as 


“B «A (8+1) (2а8-15) 
Әт 791 (B-i)(28414) (11а) 


EB A 1d, kk, "OLEI 

dg = аа + ^| рау ара (119) 
= А 

~B ~A 19 kk -1) + ilv 

YT | (5-1) (B+ p 


k. 2 
where t2= kp + ks2, p^ = z2 ‚н E (1+04К), and v = (14o-aK). 


Plane-Wave Analysis for Contraction Section 


The main stream will be regarded as compressible and inviscid 
throughout the contraction section (station B to station C of fig. 1). 
In the case of supersonic test-section flow, the term "contraction" is 
retained for convenience. As before, the turbulent field is taken to be 
incompressible and inviscid. The contraction section has its initial 
breadth and height reduced by the factors 152 and tz, respectively, 


while the velocity U(x) at station B is increased to 14U(xi) at 


station G. A cubical fluid volume element of edge D at station B will 
have been distorted into a parallelepiped of edges 14D, 12р, 15р upon 


reaching station C (fig. 5). The effect of a contraction upon a tur- 


bulent field arises principally from changes in vorticity following 
such distortion of the fluid elements passing through the contraction. 


At station B (time t=0) in figure 3, a particle at distance x 
from & corner particle of & given fluid element will &t station C 
(time t=t) be at a distance X from the corner particle. The coordi- 
nate axes are taken to move with the main stream at velocity U(x})- 
With the assumption of a weak turbulence field, the relative displace- 
ment of adjacent particles in a given fluid element due to turbulent 
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mixing is taken to be very much smaller than the displacement due to the 
contraction. The relation between x and X is then simply 


Xy = Xy (12) 


With equation (12) the continuity equation for the main stream in 
Lagrangean form provides the relation 


0111,1. = 1 (15) 


where o is the ratio of stream density at station C to stream density 
at station B. The product lolz represents the ratio of tunnel cross- 


sectional area at station C to tunnel area at station B. The parameter 
by represents the speed ratio referred to these stations. 


The equations describing the changes in vorticity following dis- 
tortion of a fluid element are, from reference 9: 


X 
af =o) a? зу 


Use of equation (12) linearizes these equations relating the upstream 
and downstream vorticities to 


E 


ш, ы отоу? (14) 


Upstream of the contraction at station В, a single plane wave being 
carried along by the main stream is designated at time t= 0 by 


POR 


i^ - Que е1Е'Х (15) 
The vorticity at station B is obtained from the curl of equation (15). 


A velocity distribution at station C compatible with equation (14) and 
satisfying continuity, equation (2), is obtained in reference 3 as 


е 


с_до iex 
q^ = 40 е 


where the wave-amplitude vector is 


JAUNES Чы Fe r (16) 
ы; dy y. s oz 
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and where the new wave-number vector K resulting from distortion of 
the fluid volume element is given by 


ME з 


Ket 2 3 
T Ir 12° 15 


(17) 


Thus both the wave-number апа wave-amplitude vectors of а plane wave ате 
altered in going through a contraction, whereas only the amplitude vec- 
tor is altered in traversing a screen. 


Equations (16) and (17) describe the effect of an arbitrary con- 
traction on a representative plane wave. For an axisymmetric contrac- 
tion defined by the condition 12 = lz (all cross sections are similar 
but not necessarily circular), equation (16) with the aid of equation (2) 
simplifies, in expanded form, to 


a,° = т —— (18a) 


Q (1-е) 
pe e B 
Ea шы 
2 ЕК + t 
~ В 
«C 1 [= в, % 5k) 
= — + ت‎ SÑQ QƏMOOA (18c) 
a5 12 a5 ek? + t? 


where € = 1,2/ не. For an axisymmetric contraction, the contraction 
parameters 17, lo, and € may be expressed in terms of the Mach num- 
bers at stations B and C as follows: 

2 

>) | 


1,2 = E ) (e) (19) 


2 
OE) 
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Spectral Tensors for Multiple-Screen-Contraction Configurations 


Equations (11) and (18) describe the effect of a screen and an 
axisymmetric contraction (152-12), respectively, upon the amplitude 
vector dy of a single plane wave typical of the assembly of waves 
representing the turbulence field (equation (1)). In the Fourier inte- 
gral Чү corresponds to -dg,, Q (K) to Q, dk,dkpdk,, and Q (K) to 
Өү dKjdK5dKz. Since at station C (fig. 1) the distortion resulting from 


the contraction transforms the wave-number vector from k to K and 
that for axisymmetry dk dkpdkz = 11122 dKjdK,üK;, equations (11) and 
(18) yield : 


iQ. Pk к, 2 
В _ A I I Bla-1)7 + цо) 
СЕЕ" t? (В-1) (28+11) (20B) 


19-Ак kz 2 
B A L 1 B(a-1)* + i(v ) 
Qz = Qz + — _ | ocn pti 2B4iu J (20c) 


2 2 
k^ + t 
2 
ai’ = 12 qi? — (20а) 
Ek] + € 
n B Qr xk; (1-6) 
ЕК + 
C p Qy ks (1-6) 
Qs = و11‎ |95 + (20f) 


єк? + t2 


If the fluid element volume T is taken to be a cube of edge D 
at station A, and hence at station B, the volume will have been dis- 
torted into а parallelepiped of edges 17D, 120, 120 at station C for 


an axisymmetric contraction. The energy spectral densities which enter 
directly into the calculation of turbulence fluctuation velocities are 
obtained from equation (3) as 
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AI 
1 2 


With the products Qa," from equations (20) formed and with the use of 
equations (21) and the continuity relations Qk. = О and Q "Kk, = 
the energy spectral densities may be written as: 


2, 2 2,2 
г (к) = a Waal г ^(x) (22a) 
шг акти / 11 | 
2 2.2), 2 
[Го2(® + EST = a T, (5) + DG) + EISE P A(x) 
4k] + pu б 
(22b) 
2 „2\2 
po „38. = k. гуу (к) (220) 
11 1 ek 


21-2 2 2,p2 
C B , |ку°(1-є)°5° - 2k (1-e) (ek) t |, в 
[5200 + ENG] = 2 = + Тыю] + j PENE ri, (8) 


(22а) 


With the use of equations (22a) and (22c), the longitudinal energy spec- 
tral density at station C for N screens in series followed by an 
axisymmetric contraction may be expressed in terms of the spectral den- 
sity at station A as 


2 2 2,2NN 


D 22 V. (ap? k. ^ev 
C 2 [Ку 4a ky +u E” 


dis conciseness, equations (228), (22b), and а шге written ав 
H.P = ABA, vB = ey. + DBA, "x = 11V1? x oe p 
Then for N screens in series, Hy = АНЕ = SE”, Vy Vi y + Dy 1? 


апа ag = hu + Pg. The lateral energy spectral nd. at 


station C for N screens in series followed by an axisymmetric con- 
traction may then be grouped as 


12-22) k 2 y AA 
C C C C 1 1 (1-e) 1 A 
[г 22 (К) +055 «| = aoe (К) + Tag (К) "rui Ak, htt? RP ER 4121262 Ea (s) es) 
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Equations (23) and (24) relate the energy spectral densities downstream 
of a multiple-screen-axisymmetric-contraction configuration to the cor- 
responding upstream spectral densities at station A. 


Results for Negligible Decay 


The solutions to be given (see appendixes B and C) will now be 
restricted to the case of isotropic upstream turbulence. The upstream 
energy spectral densities Г (к) may then be obtained from equa- 
tions (4). 8 


Turbulence velocity ratios. - As shown in appendix B, the turbulence 
velocity ratio or ratio of mean-square fluctuation velocities downstream 
of a series of М identical screens followed by an axisymmetric con- 
traction to the corresponding upstream fluctuation velocities is given 
for initially isotropic turbulence by 


C 
2 3( 
91 /N за A" sin? ө аё 


н A sint ө аб (25) 
«y A Јо (8 - cos” 9) 
C C — с 
(yn a has tata |“ ass? 0 cost om aay 
YG а G asa ый nt co g 


2 2 2 2 
2 = 1 ب„‎ Аа cos”0 + Vv“ gin^O 


4 cos^0 +u sin Ө 


A convenient approximation for equation (26) is presented later (see 
equation (39)). 


For М = 1, equation (25) for the longitudinal turbulence velocity 
ratio integrates to 


ا 
g , u PEN s ? am gsm 2 а ak tanh” E ra? tann 2‏ 
a a, a‏ 


(a2) 414 u^ E (a? -12)2 7] 


(27) 


 —  — M —‏ — کد 
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2 
12 = -Б— 
pea 
2.2 2 2 2 
Ау = а2(а2+1) + 12(1-5а2) + 8^-1)^ + (аё+1) (72-1) 


дё 


2(12-1) (n?- £7) 


A2 


Equation (26) for the lateral velocity ratio (see appendix B) 
integrates for N= 1 to 


(22), 2 252 
l _ @ ven L tanh“? i) " 1 -1 і 
a у 812 tu Ë +22 t a/ SN en 7 
(28) 
where 
: 2 
Biz РЖ (E 2-92) (2-51) s 6(Е2-12) өй а — 
2 чуг ? 
d em кары 
а _ 2. 2) 4(a*-n^) 2 


For the case of axisymmetric contraction with the screen absent 
(а2-1, K- О), equations (27) and (28) reduce, respectively, to 


(22), 0 ا‎ ) Sa^ Pre 3] 
qi? MED EU pem ше ) 


ала 


TE MUN ses [e +1) - (a? 17 (5 tanh“! 2) 
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which, in the present notation, are identical with the corresponding 
results of reference 3. Similarly, for the case of a screen and no 
contraction (а2 =» œ), the results of reference 2 are recovered in the 
form 


2t^y^ 7 


LAC _ АВ 
Eh. Bt fst? = 2+ за) - (n? 204 DIT 
1 


(2), ). Ps: ). Em га s spe 2+2) _ g]- 5 зол (nê -1) (n@-t A(z tom — i^i) 


GS) a 1st 


Punched-card equipment was used to obtain the turbulence-velocity 
ratios listed in table I. For the cases N= 2, 3, and 4, the integra- 
tions required for equations (25) and (26) were performed numerically 
by use of Simpson's rule after changing the variable of integration 
from 0 to x by applying the transformation x= cos @ Intervals 
Ax = 0.01 were used in the range О <x < 0.9; intervals Ax = 0.001 
were used in the range 0.9 < x < 1. 0. In all computations the Mach 
number Мв upstream of the “contraction was taken equal to 0.05. The 
turbulence velocity ratios listed in table I may be corrected for values 
A Mg Other un 0.05 as follows: Values of the parameters 14° , 

122 y and а? = —— for Mp = 0.05 and for the desired value of Mp 


l-e 
are obtained from equations (19). Noting that the quantities 
2 21% )y 2 v ч 2 
11 т and lo —— A depend only upon a“ and K, the values 
q 92 


of these quantities for the ac corresponding to the desired Mg are 
obtained from table I. With 142 and 15^ known for Mg = 0.05 and 
the desired Mp, the corrected velocity ratios are obtained by simple 
computation. The following empirical relations (reference 1) were 
utilized in obtaining numerical results: | 


2 
22 (EX) for K< 1 


I 


8+K 


2 _ [1.21 
a - (а) for K> 1 


For design purposes, the screen pressure-drop coefficient K may 
be estimated; according to reference 10, from the solidity ratio b, 
where b is the area of the holes in a unit area of screen, as 
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„l-b 
К x ы ^ 
For square-mesh screen with wire diameter d апа mesh designation m, « 
the solidity r&tio as defined is 
b = (1-mà)^ 


A better agreement with the screen data given in reference 1 is obtained 
from | 


bz (1-ma) 7/4 


The variation of the longitudinal and lateral root-mean-square 
velocity ratio with speed ratio 11 for a single screen (N=1) upstream 
of the contraction is plotted in figures 4(a) and 4(b), respectively, 
for selected values of the screen pressure-drop coefficient K. The 
results for K = 0, which correspond to the case of stream convergence 
or divergence in the absence of any screen, are, of course, identical 
with the results of reference 3. In general, both the longitudinal and 
the lateral fluctuation velocities downstream of the screen-contraction | 
configuration are reduced as the screen parameter K is increased. The | 
somewhat anomalous trend of the longitudinal velocity r&tios for values 
of the speed ratio less than 2 seems to reflect tbe variation of the 
auxiliary screen parameter te which approaches zero at К = 2.76, 
becomes infinitely large in the negative sense as K increases to 5.28, 
and becomes infinitely large in the positive sense as K decreases to 
5.28. 


The losses incurred through the use of damping screens are propor- 
tional to the product NKU,?, where N denotes the number of identical 
Screens in series (multiple screens) and NK is the over-all screen 
pressure-drop coefficient. The velocity ratios for a multiple-screen 
arrangement upstream of a contraction are compared on the basis of equal 
Screen losses in figure 5 for the particular case NK = 6. The advan- 
tages of using a number of screens in series to attain a given over-all 
coefficient NK are obvious. An examination of table I indicates that 
the use of multiple screens to attenuate the downstream fluctuation 
velocities becomes more effective as the over-all coefficient NK is 
increased. The screen losses can be reduced by decreasing the settling- 
chamber stream velocity Ua. Low-turbulence wind tunnels are generally 
characterized by their many damping screens and large-cross-sectional- 
area settling chambers. 


One-dimensional spectra. - In accordance with equation (6), the 
one-dimensional spectra at stations A and B are given by 
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к^ = г | г) dkodkz (29) 


E = 2 MENS акак» 


As pointed out in reference 3, а comparison of the upstream апа down- 
stream spectra on the basis of the upstream longitudinal wave number kj 
is equivalent to a comparison of the time spectra indicated by fixed 
hot-wire probes located at the corresponding stations. Defining the 


downstream spectra F“ (k) = пас such that 
со 


and 


© 
E," (ку) dk, = | E," (ку) dk, the one-dimensional spectra at 


station C are given by 


C 2 с[К1 Ez ks 
rl, ا‎ А atts Go 


Evaluation of equations (29) and (30) requires that the amplitude func- 

tion G(k) in equation (4) be specified. Compatible with the empirical 

relation for isotropic turbulence obtained in reference 8, this function 
is taken to be 


H 
G(k) = —— (31) 
(ку2+п2+62 )5 
where the constants п and Н are defined ав n= 1 1 
Мы АЁ (14) 
= 2n[ 2) . 
йе 212) 


As shown in appendix С, the one-dimensional spectra obtained from 
equations (4) and (31) may be expressed in terms of a dimensionless wave 
number k, /n as incorporated in the following parameters: 


l! 


8 l + k/n? 


|! 
m 
— 
= 
Do 
+ 
TN 
= 
№ 


г = 
(52) 


0% 
| 
U 
— 
г 
Do 
+ 
Š, 
A 
Yo 
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Thus the upstream one-dimensional spectra for this special case of 


isotropic turbulence are, in dimensionless forn, ` 
A 
F. (k/n) š 
ee (33) 
A 8 
Fi (0) 
A 


Fo“ (0) 52 


Also, the longitudinal one-dimensional spectrum downstream of a single- 
screen-axisymmetric-contraction configuration may be written (see | 
appendix C) in dimensionless form as 


C à 
1 2v s+h 4(s-1) 
1 2 5 
F4A( 0) a е 8 e u2 E 


2 2| 2,2 , 2 2 
c = ES сш " gh Р e es " Cm) u |а. E m je у 
h f S v a^(u^-4) - u 


2 2 
С. = 3sg , g(s+2) + 2s _ „ _ ue (sth) (eh) 
2 Wa 2 “ше E NE MEI E DLL ale2(u2-a) - y2] 
с. = (83-2) (£-8) (u*-4)2a# 
ETICA Е ка 


The corresponding lateral one-dimensional spectrum is 


C 

Po (kı/n)  2(s.3) sth 4(в-1) 
A - 2 |E] + F, 108, — 5 108. a 

Fo (0) 12 8 


where 
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= a^ (58-2) M" (v*-a*u*) (28-2) ñ 222 [g(h-f) g - 1] 
- 2s^(s-1) 2u°Pes a“fh fh 28 
o f(grh)-g(h-f) , ula" (v?-4a7) - v^ 2h(8-£) -fg 
Zz, | 2(sth) 2-2 + 2.272 22 af + 2 
2afuĉh f^h yn“ a? (2-4) =p f^h 
2.2 (sth) 202-402) - v2] on + 558 , 4(а2-1) (02-а2р2)һ 
2 а2ъ4 2012-4) - p2] a [82(u2-4) - uZ][a2(v2-4a2) - vZ] 
— 4(в-1) (v^ -a 4,7) Г. 4 ° 
Ez = 1 + — ——35 
H a^(u^-4) - u 


The one-dimensional spectra given by equations (33) to (36) are 
applicable when the amplitude function G(k) . has the particular form of 
equation (31). Although these spectra are not expected to be valid for 
the very high wave numbers because of the neglect of viscosity, various 
experiments on isotropic turbulence have indicated that equations (33) 
and (34) provide a very good approximation to that portion of the actual 
isotropic spectrum containing the largest part of the turbulent energy. 
Equations (35) and (36) should furnish a similar approximation for 
axisymmetric turbulence. The restrictions given for equations (33) to 
(36) do not apply to the expressions for turbulence velocity ratios, 
equations (25) and (26), for which there is no need to particularize the 
spectrum amplitude function G(k). 


The downstream longitudinal and lateral one-dimensional spectra, 
equations (35) and (36), are compared with the corresponding upstream 
isotropic spectra, equations (33) and (34), in figures 6(a) and 6(b), 
respectively, for the following typical case: Mp = 0.05, Мс = 2.0, 


К = 2, М = 1. The case К = 0, as obtained in reference 3, has also 
been included for comparison. Тһе scaling factors indicated by equa- 
tions (B5) and (B6) of appendix B have been incorporated in the down- 
stream spectral ordinates so that the zero-wave-number intercept gives 
the turbulence scale ratio (appendix B). 


The distortion in shape of the longitudinal spectrum noted in ref- 
erence 3 as a consequence of the stream convergence is accentuated 
(fig. 6(a)) by the presence of a damping screen upstream of the contrac- 
tion. This distortion is accompanied by a reduction in the ordinate 


— C — A —2\С —5\В 
values by the factors (5,2) m ) amd (ч.2) 3 for K = 2 
and К = О, respectively. The downstream lateral spectrum ordinates 
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——\С//—=\А 
(fig. 6(b)) are increased by the factors (92) (92) and 
do^ / 422) for К= 2 and КЕ = О, respectively. The distortion in 
shape is relatively slight compared with the distortion noted Рог the 
longitudinal spectrum. 
As may be seen from equations (33) and (34) for the upstream iso- 


tropic spectra, the longitudinal and lateral spectral ordinates have 
maximum values at к/а = О and kj/n = 1/A/3, respectively. The 


situation is reversed for the downstream spectra. Here the lateral 
spectral ordinates have maximum values at kj/n = О and the longitu- 


dinal spectral ordinates at k,/n% 1.4. Occurrence of a peak in the 


spectrum curve at some wave number other than zero is an indication 
that the correlation coefficient may take on negative values. 


Scale ratios and correlation coefficients. - For the scales of tur- 
bulence defined by equation (10), the longitudinal and lateral turbulence 
scale ratios (ratios of downstream to corresponding upstream scales) for 
a screen-contraction configuration are obtained in appendix B as 


zr 


C N =Ë 
Cw гу (a? 
o | tz ее 


43 
C maak Ol aan 
ae an p a (38) 
2 


The scale ratios obtained from equations (37) and (38) which do not 
require that the amplitude function G(k) of equation (31) be specified 
are listed in table I for the case of isotropic turbulence at station A. 
Typical results are plotted in figure 7. 


The lateral scale ratio (see fig. 7(a)) approaches a constant value 
of approximately 4/3 for values of the speed ratio 14, greater than 3. 
Measurements of the lateral correlation curve at a speed ratio near 
unity which are reported in reference ll indicate that the lateral scale 
is substantially unchanged by damping screens. This is in qualitative 
agreement with the present result which indicates that for 11 slightly 


greater than unity the downstream lateral scale will not exceed the 


im Р-1 -kj/n 
or example, when E, (ка) = Ку e ; the correlation coeffi- 


cient is obtained by using equation (7) as 
1\2 -P/2 
(2) + гу Г (Р) сов (Р фап-1 nr, ) where T designates the gamma 


function. For Р = 1, Ry(Fڊ)‎ is always positive; for P> 1, В. (гл) 
will take on negative values for particular values of rj. 
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corresponding upstream scale by more than about 20 percent. Taking the 
lateral scale ratio equal to 4/3 leads to the following convenient 
approximation for the lateral turbulence velocity ratio from 


equation (38): 
=)" 
" үн E s eN 
= eo E, 


For a given value of the screen pressure-drop coefficient МК, the 
longitudinal scale ratio (see fig. 7(b)) decreases with increasing speed 
ratio 14 to а minimum value at 14 = 27.4 (corresponding to Mp = 0.05, 
Мс = 4/3) where the contraction parameter a? has its minimum value. 
As shown in table I, the longitudinal scale ratio attains a zero value 
when the screen parameter v^ = O (К #2.76). This and the occurrence 
of maximums in the downstream longitudinal spectrum curves at nonzero 
wave numbers suggest that the downstream longitudinal correlation coef- 
ficients are negative for extensive ranges of the separation distance 
r4. Under these conditions interpretation of the conventionally defined 
scales as lengths characteristic of the average size of the turbulence 
eddies is open to question, and consideration of the correlation coeffi- 
cient curves is advisable. 


(59) 


The correlation coefficients at station А Рог isotropic turbulence 
with the spectrum amplitude function G(k) given by equation (31) are 
obtained from equation (7) as 


Ry =. 3P (40) 


r-n 
A 9 -r,n 


The contour integrations used to obtain equations (40) and (41) are not 
valid when гу = 0; hence the microscales Ay are evaluated from the 


integral relation of equation (9). Such evaluations indicate that 
a = Ns O, which is to be expected in view of the neglect of vis- 


cosity effects in the analysis. The longitudinal correlation coefficient 
curve of equation (40) is plotted in figure 8(a) and is always positive; 
the lateral correlation coefficient curve of equation (41) plotted in 


figure 8(b) reaches its zero value at гуп = 2 (ру = 214 ) and its 
minimum value at ryn= 3 (ху = arar 


The downstream correlation coefficient curves (at station C) have 
been obtained numerically for the case N = 1 from the following 
rearrangement of equations (7): 


A — — —— — — — ---- 
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с \ یر‎ 
R4" (гүп) = E | ES e cos kyr, d (2) (42) 
Ro” (r4n) = 1 iis 0 E) (2) cos kyr, à (3) (43) 


Fo Ato) m ape) 


In evaluating equations (42) and (43), values of the integrand were 
obtained for k,/n ranging from O to 50; and for k,/n greater than 
ry (3) ces 
50, in view of the asymptotic behavior of the functions 
F (o) qe 
cos kr] 
the integrand was approximated as u— 275 Typical downstream a 
(к/а) 
tudinal and lateral correlation coefficient curves (for the case 
МВ = 0.05, М = 2.00, K = 2, N= 1) are also plotted in figures 8(а) 
and 8(b), respectively, to indicate the changes resulting from passage 
of initially isotropic turbulence through a given screen and contraction. 
Although the downstream lateral correlation coefficient is shown in fig- 
ure 8(b) to reach slightly negative values, it is believed that these 
are the result of unavoidable "round-off" errors in. computation of the 
Fourier transforms and that the coefficient is actually always positive, 
consistent with the corresponding spectrum curve of figure 6(b), which 
has its maximum value at zero wave number. 


The correlation between simultaneous fluctuation velocities at two 
points a distance rı apart will decrease more rapidly with increasing 
values of r4 when the eddies comprising the turbulent field are small 


than when the eddies are large. Figure 8(a) thus indicates that the 
longitudinal scale of an initially isotropic field of turbulence is 
decreased by passage through the particular screen-contraction config- 
uration chosen. Figure 8(b) indicates that the corresponding lateral 
scale is increased. 


in view of the negative values attained by the downstream longi- 
tudinal correlation coefficient, no physical meaning can be assigned to 
the longitudinal scale ratio defined in the conventional manner by 
equation (10). For example, the longitudinal scale ratio reaches a zero 
value even though the longitudinal turbulence velocity ratios are finite 
when the screen pressure-drop coefficient NK has the value 2.76. The 
negative values attained by the upstream lateral correlation coefficient 
do not present a similar anomaly because of the relation between the 
longitudinal and lateral scales in the case of isotropic turbulence, 


namely, L^ = 215^. 
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The difficulty is removed if an effective longitudinal scale Гл! 
is defined as the positive area under the corresponding correlation 
curve. Effective longitudinal scale ratios are plotted in figure 9 and 
show a qualitative similarity with the conventional ratios shown in fig- 
ure 7(b). For a given value of the screen pressure-drop coefficient 
NK, the effective scale ratio decreases with increasing speed ratio 14 
to a minimum value at 1] = 27.4 for which the contraction parameter 


а? has its minimum value. For a given contraction the effective scale 
ratio reaches its minimum value when МК ZX 2.76. 


ESTIMATION OF DECAY EFFECTS 


In view of the assumptions of inviscid flow and small fluctuation 
velocities relative to the main stream, the preceding analysis is 
strictly applicable only in the &bsence of the turbulent decay processes 
(viscous dissipation and turbulent mixing). For many wind-tunnel con- 
figurations, effects of decay upon turbulence are of the same order of 
magnitude as the screen-contraction effects.  Carrection of the theo- 
retical turbulence velocity ratios may therefore prove necessary for 
practical applications of the theory. 


Selection of the appropriate decay correction presents certain 
difficulties inasmuch as there is a lack of experimental investigations 
of axisymmetric turbulence decay. Some guidance may be obtained from 
the theoretical studies of Batchelor (reference 4) and Chandrasekhar 
(reference 12) on axisymmetric turbulence. The time rates of change of 
the mean-square velocity components are, in the notation of reference 4: 


a (u2) = -am + 20(- 10a - 2b - 2C - 14d) 
2 (u2) = em, + 2v(- lOa + b - 3d) 


In these equations and in equations (44) and (45), the symbol v rep- 
resents the kinematic viscosity coefficient. The corresponding expres- 
sion for the mean-square resultant velocity is 


E: пу + 2002) = - 2v(30& + 2c + 20d) (44) 


For isotropic turbulence, с = а = O and equation (44) becomes 


2:57) de Qu) = азон 6 
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The velocity components wy“ and Up" of reference 4 are identical with : 


ay? and ъё in the present notation. The quantities a, b, c, and d 
in appropriate groupings represent the coefficients in the series expan- 


sions in r, for the longitudinal and lateral velocity correlation 
coefficients. The quantity по depends on the two-point velocity- 


pressure correlation which tends to zero as isotropy is approached. For 
the decay of isotropic turbulence in a constant-area channel during the 
initial period wherein both inertia and viscous forces are of importance, 
equation (45) leads to the semiempirical relation (reference 13) 


PE | un t? 
922 _ + [бы mJ (46) 


cy F MOL 


where аё апа 952 represent the mean-square velocity components at 


any station downstream of the reference station A and t(11) represents 
the appropriate decay time. 


The absence of the velocity-pressure correlation term in both z 
equations (44) and (45) suggests that, provided the quantity (2c + 20d) 
is much smaller than the quantity 30a, equation (46) may yield a sat- 
isfactory approximation for the decay of the mean-square resultant tur- 
bulent velocity in axisymmetric turbulence. The data of references 1 
and 14 tend to support such an approximation. The assumption that the 
effects of the screen-contraction combination and the decay upon the 
turbulent velocity ratios proceed independently (see reference 2) leads 
to the relation 


|| (a? )s (ey N _ E. (че) ) 
ОО R^ A 


where the subscript sc refers to the turbulence velocity ratios 
obtained in the absence of decay, computed from equations (25) and (26) 
and listed in table I, and the subscript scd implies that the effects 
of initial period decay have been included. In computing J from equa- 
tion (46), the decay time t(17) is taken as the time required for a 
particle at local main-stream velocity to pass through the screens and 
the contraction starting from station A. This implies that the contrac- 
tion affects only the decay time. Some question exists as to the appli- 
cability of equation (46) and hence equation (47) for damping screens 
in which the wire diameters are usually very small. 


(47) 
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A comparison of the theoretical mean-square resultant turbulence 
velocity ratio corrected for decay by the use of equation (47) with the 
experimental ratios obtained from reference 1 is shown in figure 10. 

The mean-square resultant velocity in the absence of decay for the case 
of a single-screen-contraction configuration (N-1) is also included to 
show the magnitude of the correction involved for the configuration of 
reference 1. The following data were used in applying the decay cor- 


rection: UA = 62-8 feet per second, КЕ BA /2 = 0.15 foot per sec- 


ond; and L^ = 0.05 foot (estimated). The screen pressure-drop coef- 
ficients were corrected as suggested in reference 14 


„Ае 
е 2 асл 


where К. designates the screen pressure-drop coefficient measured at а 
given speed Up. Although the single experimental points obtained for 


each multiscreen arrangement do not check the decay correction as well 
as do those for the single-screen arrangement, the limited data do not 
warrant any refinement of the correction method for multiscreen- 
contraction configurations. 


ln order to obtain the resolution of the resultant turbulence veloc- 
ity ratio into longitudinal and lateral components, some knowledge of 
the velocity-pressure correlation is required. As shown in reference 4 
the effect of this correlation as represented by the term m is to 
transfer energy from the larger to the smaller of the velocity components, 
thus providing a drive towards isotropy. As shown in table I, the lon- 
gitudinal component will, in general, be much smaller than the lateral 
component so that adjustment of the longitudinal component is more crit- 
ical than adjustment of the lateral component. The magnitude of the 
longitudinal component is governed by two opposing effects.  Turbulent 
decay processes reduce this component; the drive towards isotropy tends 
to increase it. TIn the absence of any quantitative knowledge concerning 
the velocity-pressure term шо, the simplest assumption to be made is 
that the longitudinal turbulence velocity ratio may be corrected for 
decay and isotropy drive by taking &n average of the values for zero 
decay and isotropic decay or 


El al Ey 
(e2) -— É; "ON 


Consistent values of the lateral turbulence velocity ratio are then 
obtained from the longitudinal velocity ratio of equation (48) and the 
resultant velocity ratio of equation (47). 


(48) 


- ا — —— — ——— — سس‎ — — л 
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The comparison shown in figure 11 provides some estimate as to the 
agreement that might be expected between the predicted turbulence- 
velocity-ratio components (corrected for decay) and the experimental 
values. The agreement shown is considered satisfactory for most engi- 
neering applications. The theoretical velocity ratios obtained in the 
absence of decay are included for the case N= 1 to indicate the 
magnitude of the correction. 


The turbulence scales are also affected by the turbulence decay 
process, tending to increase as the decay time is increased. Under the 
action of the viscous forces the smallest eddies are dissipated so that 
the average eddy size (scale) would be expected to increase. For 
isotropic turbulence, the change in scale during the initial period 
analogous to the relation given for the fluctuation velocity, equa- 


tion (46), is 
2 
eon 
Lo 


Presumably the effect of decay upon the scales of turbulence could 
thus be obtained by a procedure similar to the one suggested for the 
fluctuation velocities. In the absence of any experimental data such 
development does not appear warranted. 


CONCLUDING REMARKS 


The present analysis treats, in the absence of turbulent decay 
processes, the combined effect of'a series of identical damping screens 
followed by a stream convergence (or divergence) upon the mean-square 
fluctuation velocities, scales, correlation coefficients, and one- 
dimensional spectra of a field of turbulence convected by a main stream. 
Numerical results are presented for the case of upstream isotropic 
turbulence. 


The limited experimental data available confirm at least qualita- 
tively some of the theoretical results obtained such as the distortion 
of an initially isotropic field of turbulence by the damping screens and 
stream convergence into a field axisymmetric about the main-stream 
direction with the lateral components of the resultant fluctuation veloc- 
ity larger in magnitude than the longitudinal component, and the relative 
insensitivity of the lateral scale of turbulence to damping-screen and 
stream-convergence effects. The beneficial effects of using several 
screens in series to attain a given over-all screen pressure-drop coef- 
ficient in attenuating the fluctuation velocities are also substantiated. 
This attenuation is accentuated as the screen coefficient NK is 
increased. 
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The theory predicts certain marked changes in the ordinates of the 
downstream one-dimensional spectra and, in the case of the longitudinal 
spectra, a noticeable distortion of shape which should be confirmable by 
experiment. The longitudinal downstream correlation coefficients attain 
negative values over a large range of the separation distance т]. 

Under these conditions, the scales of turbulence ав conventionally 
defined cannot be regarded as representative of the average eddy size. 
Accordingly, the longitudinal scales have been redefined. The effect 
of the damping screens and stream convergence is to decrease the lon- 
gitudinal scale and to increase the lateral scale. 


An approximate method of correcting the predicted turbulent fluc- 
tuation velocities for the effects of turbulent decay is presented. 
Tabulations of the fluctuation velocities over a wide range of conditions 
are provided for convenience in engineering applications. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, October 28, 1952 
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APPENDIX A 


SYMBOLS 


The following symbols are used in this report: 


Al ‚2 


mae 


parameter groupings defined after equation (27) 
auxiliary contraction parameter, аё = 1/1 - € 


parameter groupings defined after equation (28) 


solidity ratio of damping screen 


parameter groupings defined after equation (35) 


edge lengths of volume within which the turbulence field is 
defined | 


wire diameter of damping screen 


parameter groupings defined after equation (36) 

Fi, Fo, or Fz 

one-dimensional longitudinal spectral density (see equation (6)) 
one-dimensional lateral spectral densities (see equation (6)) 


auxiliary wave-number parameter, f = 8/12 + 4/12 


amplitude function in isotropic spectrum tensor (see equa- 
tions (4) and (31)) 


auxiliary wave-number parameter, g = s/ t^ + 402/2 


constant appearing in amplitude function of special isotropic 


spectrum tensor, H = = (q2)^ (see equation (31)) 
x 
_ 1 - a? - gl 
auxiliary wave-number parameter, h = XXE > SEES 
a 


т 
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Р 


| 


turbulence decay factor (see equation (46)) 


Ap 
screen pressure-drop coefficient, Ке T7 o 
2 


amplitude of vector К: К = куё + kọ + kz? 
E = к, Ко, ог К; wave-number vector 


і; Lo; or Ls 


longitudinal scale of turbulence (see equation (10)) 


effective. longitudinal scale of turbulence 


lateral scales of turbulence (see equation (10)) 


stream velocity at station C divided by stream velocity at 
station B (see equation (19)) 


stream breadth at station C divided by stream breadth at sta- 
tion B (see equation (19)) 


stream height at station C divided by stream height at sta- 
tion B 


stream Mach number at station B, C 


mesh designation of damping screen (reciprocal of center-to- 
center distance between neighboring wires) 


number of screens in series (cascade) 


constant appearing in amplitude function of special isotropic 
1 


ea 


spectral tensor, n = 


у^ 


constant 


static pressure 


Q. = Qi; Q. ; Or Qz ; wave-amplitude vector 


ge ee 


30 


У nu 
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dy = dj» 92, OT 95; turbulence-velocity-fluctuation vector 


correlation coefficient (see equation (7)) 


correlation tensor, Rys(r) = q (x) ag (xt) 


гү = Гү, Гр, OF F3; separation vector 


wave-number parameter, в = k,°/y* + 1 
time 
decay time 


main-stream velocity 


longitudinal component of combined turbulent velocity fluctua- 
tions and potential-flow velocities induced by screen 


longitudinal root-mean-square turbulence velocity ratio (used 
in table I) 


lateral root-mean-square turbulence velocity ratio (used in 
table I) 


lateral components of combined turbulent velocity fluctuations 
and potential-flow velocities induced by screen 


Xy = ху, X2, OT Xzj position vector 


ES 
2 2 2) 
E (к, + Ez 


three-dimensional spectral tensor 


Ae ада” 8 + v? in^ 0 
4 cos? 0 + ue sin? Ө 
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5.5 Kronecker delta; 5,5 = 1 for ry = Š and &5=0 for YF #6 
€ axisymmetric contraction parameter, є = 1,2/1,2 
à ue 
n auxiliary sereen parameter, nE = 5 
U“ = 4 
6 polar angle (see appendix B) 
K amplitude of vector K; K^ = ку + Ko“ + Kx" 
= K Ks: = 

K i К-› 22 OF ^: wave number vector at station C 
A Е REC UT (see equation (22a)) 

4k, + реб 
Ay Ay? A2» or ^z 
№ longitudinal microscale of turbulence (see equation (9)) 
№ ,3 lateral microscales of turbulence (see equation (9)) 
ш auxiliary screen parameter, р = 1 + а + K 
v auxiliary screen parameter, v =1 + a - aK 

v 
t^ auxiliary screen parameter, £2 = 5 Е 
vr ~ da 

p stream density 

(02 -а2р2) ку“ 
> — O (see equation (22b)) 

4k] + u^t 
G main-stream density at station C divided by main-stream density 


at station B 


T volume 


ee 
ar —— — — —— — -  — T  — — Rn  — A ت س‎ 
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Ф azimuth angle (see аррепдіх В) 
Xy = Хі, Ха, ог Xs; position vector (see equation (12)) 
Ya angle to screen normal of flow incidence upstream of screen 
Yo angle to screen normal of flow emergence downstream of screen 
2 с да ti Я (see equation (22d)) 
(ek, 244%) 
e 4, = 04, 02, ОГ 03; vorticity vector 
Superscripts: 
A station upstream of screens 
B station downstream of screens and upstream of contraction 
C station downstream of contraction 
* complex conjugate 
Subscripts: 
A station upstream of screens 
B station downstream of screens and upstream of contraction 
C station downstream of contraction 
N number of like screens in series 
sc only effects of screens and contraction present 
scd effects of screen and contraction corrected for initial period 
of decay 
1 longitudinal component 


lateral components 
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TURBULENCE VELOCITY AND SCALE RATTOS 
Velocity ratios. - Using spherical polar coordinates 


Ку =k cos Ө 
ko =k sin Ө cos q 
kz =k sin Ө sin @ 
кё = je, Š + k^ + к^ 
equations (25) апа (24) may be put in T form 


2,2 
lo k G(k) AN sin? ө 


TC ср эы шо NN Bl 
E 11 ЈЕ 11 (€ cos? Ө + sin? 0)2 (B1) 


2,2 N-1 ¿2 2 2 


(4 cos? Ө + u2 sin? 0)(c cos? Ө + sin? ө)? 


442 cos? 9 + v? sin? 8 
4 cos? 0 + ue sin? Ө 


where AS 


The downstream mean-square fluctuation velocities are given by 


a?) = Irse aK aka 


analogous to equation (5). Inasmuch as the function G(k) appears in 
the expressions for the energy spectral densities, the variable of inte- 
gration will be changed from K to k so that 


8 
— 1 
ay? = > J Uy (g) dkjdkedkz. Noting that 


112^ 5 


dk, dkodkz = k^sin Ө дө аф dk, the downstream mean-square velocity compo- 
nents of the turbulent field are obtained from equations (Bl) and (B2) as 
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and, inasmuch as the downstream turbulence will be axisymmetric when 
the upstream turbulence is isotropic, 


a, e (2), n аё (e. + Sot I PY dk Jo ao AH-l gin? Ө cos? Ө 10 
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(4 cos? Ө + p2 sin? 0) (є cos? Ө + sin? 0)? 


The mean-square velocity components of the upstream isotropic turbulence 
are obtained by using equation (4) as 


"ol = "3l Е (a) E de k^G(k) ак E" аф o sinje do = = Jo xem) Е 


The turbulence velocity ratio or ratio of mean-square fluctuation 
velocities downstream of a series of N identical screens followed by 
an axisymmetric contraction to the corresponding upstream fluctuation 
velocities is then given for the longitudinal and lateral components, 
respectively, by 


(2), NU ИМ _ — AN gin? Ө 40 (вз) 
( "OMNE (a2 - cos? Ө)? 
91 
ae 
(=), Mc en. MCI bs s 3(a 1922 -x* y) АН-1 sin? 8 cos? Ө 40 (B4) 
Go (z2)^ Mb Blo (4 cos^644* sin” ө) (a^ -cos? p) 
where al = LI The new contraction parameter al is introduced for 


convenience in subsequent calculations. It may be noted that the veloc- 

ity r&tios are independent of the amplitude function G(k) which appears 
in both the isotropic and axisymmetric spectral tensors. Equations (B3) 

and (B4) appear in the text as equations (25) and (26), respectively. 
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Turbulence scale ratios. 


- The turbulence scales may be obtained 
from the energy spectral densities as indicated by equations (6) and 
(10). Compatible with the formulation 


(2,2) Y а, (к) dk, dk, ак, 


the longitudinal scale at station C is 


or, applying equation (24), 


(La), іб B A(x) | dkodk 
toa (a2) Né J 11 во i 
ta - | 


The longitudinal scale at station А is 


(1^ = yd) BO] со TAs 
91 / -œ 


The ratio of longitudinal scale downstream of a series of N identical 


Screens followed by an axisymmetric contraction to the corresponding 
upstream scale or longitudinal scale ratio is thus 
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The corresponding ratio for the lateral scales is obtained in a similar 
manner as 


21-1 
(ъа) _ вео > N (9,2 " 
pups = — = (a 122 (B6) 
(La) sk (ZY ). "o 
[F2(0) | OS 
2 


These relations for the scale ratios do not require that the upstream 
turbulence be isotropic. Equations (B5) and (B6) appear in the text as 
equations (37) and (38), respectively. 
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ONE-DIMENSIONAL SPECTRA 


With the use of equations (4) and (31), equations (29) can be written 
со 
ЛА ES at 
5 
(k 2 + n? + t?) 
"s 
(ak + €2)€ at 


pA =F," = oTH —— 


2 
О 


Integration yields, after use of equation (32), 


E 
Fy = = (c1) 
n 5 
ка = т = 03820) (а) 
218 


Equations (33) and (34) follow upon dividing equations (C1) and (C2) by 
A A 

F and (Е respectively. 

1 ) /n=0 ( 2 Jk /n=o’ 


With use of equations (4) and (B1) and (B2) of appendix B, equa- 
tions (30) can be written 


: "T (40%) 2422)" (= 22) Bat 
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(407%, 2402 ES) BUT (10,24 2) ^y, 2tŠat 


a ateo. 2a(a2-1)2(2-a2u2) | (c4) 
va ) no үг Ju- у c | (262202) (k 2.2442)" oP (2st?) - в? | 


For the case of a single-screen-axisymmetric-contraction configuration, 
integration of equations (C3) and (C4) yields equations (55) and (36) of 
the text, respectively. For the case N= 1, the quantity (F2 1 
designates the one-dimensional lateral spectrum downstream of an axisym- 
metric contraction in the absence of damping screens and may be obtained 
from reference 5. 
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Figure 1. - Configuration treated in analysis. 
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Figure 2. - Action of damping screen on components of combined turbulent and 


induced velocities at screen. " 
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Figure 3. - Typical fluid-element- and plane-wave distortions resulting from 
stream convergence. 
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Root-—mean-square longitudinal turbulence velocity ratio, 
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Speed ratio, l 
(a) Longitudinal turbulence velooity. 


Figure 4. - Variation of root-mean-square turbulence velocity ratio with speed ratio (Mp of 0,05) and soreen prassure-drop ooef- 


ficient K in absence of turbulence decay for single-soreen-axrisyomeatric-contraction configurations with upstream isotropic 
turbulence. 
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(b) Lateral turbulence velooity. 
Figure 4, ~ Conoluded, Variation of root-mean-aquars turbulenoe velooity ratio with speed ratio (Mg of 0.05) and screen pressure- 


drop рое ова К in absence of turbulence decay for single-sorsen-axiaymmetrio-contraction configurations with upstream 160- 
tropic turbulen 
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Figura 5. - Effect of multiple screens Н and speed ratio (Mp of 0.05) on root-mean-square turbulence velooity ratio in absence 


of turbulence déoay for soraen-axiayumetric-contraotion configurations with upstream isotropic turbulence and constant soreen 
losses. Over-all scraen pressure-drop coefficient, NK, 6. 
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Speed ratio, 1; 
(b) Lateral turbulence velocity. 
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Figure 5. - Oonoluded. Effect of multiple soreens Н and speed ratio (Mg of 0.05) on root-mean-square turbulenos velocity ratio in 


absonoe of turbulence deoay for soroen-axisymmetrio-oontraotion configurations with upstream isotropic turbulenes and constant 
soreen losses. Over-all вогееп pressure-drop coefficient, НК, 6. 
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(a) Longitudinal spectra. 


Figure 6. 
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- Comparison of one-dimensional spectra in absence of turbulence decay for contraction 


and for single-screen-contraction configurations for upstream isotropic ‘turbulence having 


amplitude function G(k) = H(k? + n2)-3, 


Mg, 0.05; 


Ma» 2. 00; LE 29. B22. 
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(b) Lateral spectra. 
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Figure 6. - Concluded. Comparison of one-dimensional spectra in absence of turbulence decay for contraction 
and for single-screen-contraction configurations for upstream isotropic turbulence having amplitude func- 


tion G(k) = H(k^ + п2)-5. мы, 0.05; Mp, 2.00; 1, 29.822. 
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(a) Lateral scale ratio. 


Lateral scale ratio, (12)5/(22)* 
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Figure 7. - Variation of scale ratio with speed ratio (Mz of 0.05) and screen pressure-drop coef- 


ficient K in absence of turbulence decay for single-screen-axisymmetric-contraction configura- 
tions with upstream isotropic turbulence. 
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(b) Longitudinal scale ratio. Longitudinal scale ratio equals zero for К» 2.76. 


- Concluded. Variation of scale ratio with speed ratio (Мв оГ 0.05) апа screen pressure- 
drop coefficient К in absence of turbulence decay for single-screen-axisymmetric-contraction 


configurations with upstream isotropic turbulence. 
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(a) Longitudinal correlation coefficients. 


Figure 8. - Comparison of correlation coefficients in absence of decay for a 
screen-contraction configuration (Mp = 0.05, М = 2.00, K = 2, Ñ = 1) with 


upstream isotropic turbulence having amplitude function G(k) = н(к2 + n2)73, 
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(b) Lateral correlation coefficients. 


Figure B. - Concluded. Comparison of correlation coefficients in absence of 
decay for a screen-contraction configuration (Mg = 0.05, М; = 2.00, K = 2, 
N = 1) with upstream isotropic turbulence having amplitude function 
G(k) = H(k® + п2)-5. 
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Figure 9. - Variation of effective longitudinal scale ratio in absence of turbulence decay with speed ratio (Mp of 0.08) and screen 


pressure-drop ooafficient K for single-soraen-axisymmetrio-contraation oonfiguration with upstream isotropic turbulence having 
amplitude function @(k) = H(k? + n2)-5, 
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 Mean-square resultant turbulence velocity ratio, | 
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Over-all soreen ргеввиге-@гор coefficient, NK 


Figure 10. - Comparison of theoretical mean-square resultant turbulence velocity ratios corrected for decay with sxperimant of 
reference 1. Speed ratio 1], 6.7 (Mg = 0.06, Ho = 0.34); М screens in seriesj upstream isotropic turbulence; scale 


Le^, 0.05 foot (estimated). 
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Figure ll. - Comparison of theoretical root-mean-square longitudinal and lateral turbulence velocity ratios corrected for decay with 
шеш of reference 1. Apeed ratio 1, 5.7 (Mg = 0.05, Mç = 0.54); Н soreens in series; upstream isotropic turbulence; 


scale r A , 0.05 foot (estimated). 
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